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A self-consistent renormalization (RG) scheme has been apphed to nonhelical magnetohydrody- 
namic turbulence with zero cross helicity. Kolmogorov's 5/3 powerlaw has been shown to be a 
consistent solution for d > dc ~ 2.2. For Kolmogorov's solution, both renormalized viscosity and 
resistivity are positive for the whole range of parameters. Various cascade rate and Kolmogorov's 
constant for MHD turbulence have been calculated by solving the flux equation to the first order 
• in perturbation series. We find that the magnetic energy cascades forward. The Kolmogorov's con- 

I slant for d = 3 does not vary significantly with ta and is found to be close to the constant for fiuid 

04 . turbulence. 

Oh! 

cn 
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' The statistical theory of magnetohydrodynamic (MHD) turbulence is one of the important problems of current 
research. The quantities of interests in this area are energy spectrum, cascade rates, intermittency exponents etc. In 
(— I this letter we analytically compute the renormalized-viscosity, renormalized-resistivity, and cascade rates using the 
field-theoretic techniques. 

The incompressible MHD equation in Fourier space is given by 

CN ; {-iuj + iyk'^)u^{k) = -^P^„^{k) J dp[uj{p)u,nik-p)-b^ip)b,^^^ (1) 

^ : (-^^ + Afc^) 6.(fc) - -^P,^M f dp[u,{p)b^{k - p)] (2) 

O 

m . 

O . k^h{k) = (4) 



hu.ik) = (3) 



where u and b are the velocity and magnetic field fluctuations respectively, v and A are the viscosity and the resistivity 



respectively, and d is the space dimension. Also, 



P^^M = fc,P,™(k)-t-fc„P,,(k); (5) 

>: P^m{k) = 5™-^; (6) 

■ k^ 

; PtrmO^) = kjS,r,r - kmS^f, (7) 

; k^(k,uj); dp = dpduj/{2Trf+\ (8) 

The energy spectra for MHD, E^{k) and E^{k), are still under debate. Kraichnan [Q and Irosnikov ||^ first 
gave phenomenology of steady-state, homogeneous, and isotropic MHD turbulence, and proposed that the spectra 
is proportional to k~^/'^. Later Marsch Matthaeus and Zhou [Q, and Zhou and Matthaeus ^ proposed an 
alternate phenomenology in which the energy spectra are proportional to k~^^^, similar to Kolmogorov's spectrum 
for fluid turbulence. Current numerical 1^-^ and theoretical [9|-|ll| work support Kolmogorov-like phenomenology 
for MHD turbulence. In the present paper we show that Kolmogorov's spectrum (cx k^^/'^) is a consistent solution of 
renormalization group (RG) equation of MHD turbulence. 

Forster et al., DeDominicis and Martin, Fournier and Frisch, Yakhot and Orszag |l^ applied RG technique to fluid 
turbulence. They considered external forcing and calculated renormalized parameters: viscosity, noise coefficient, and 
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vertex. McComb and his coworkers instead applied a self-consistent RG procedure; here the energy spectrum 
was assumed to be Kolmogorov's powerlaw, and the renormalized viscosity was computed iteratively. For MHD 
turbulence, Fournier et al., Camargo and Tasso, and Liang and Diamond employed RG technique in the similar 
lines as Forster et al. [|l2|. In this letter we will apply McComb's self-consistent technique to MHD turbulence. 
Earlier Verma j|] had done a self-consistent calculation and showed that the mean magnetic field gets renormalized, 
and the Kolmogorov's powerlaw is a consistent solution of MHD RG equation. In this letter we will carry out the 
renormalization of viscosity and resistivity. 

For simplicity of the calculation we assume that the mean magnetic field is absent. This allows us to assume the 
turbulence to be isotropic to a reasonable approximation. In presence of mean magnetic field, turbulence become 
anisotropic; this issue has been studied by Sridhar and Goldreich and Goldreich and Sridhar In addition to 
the above assumption, we also take cross helicity (2u • b), magnetic helicity (a • b), and kinetic helicity (u • cu) to be 
zero, where a is magnetic vector potential, and u is the vorticity. 

In our RG procedure the wavenumber range {kN,ko) is divided logarithmically into N shells. We carry out the 
elimination of the first shell fc> = {ki,ko) and obtain the modified MHD equation for k"^ = {k^jki). This process 
is continued for all the shells. The shell elimination is performed by ensemble averaging over k^ modes We 
assume that uf{k), and bf{k) have gaussian distributions with zero mean, while uf{k) and hf{k) are unaffected by 
the averaging process. In addition we take 

{u>{p)u>{q)) = P,,(P)C""(P)'5(P + q) (9) 
{b>{p)b>{q)) = P.MC'\mP + q) (10) 

Let us denote i^(„) and A(„) as the viscosity and resistivity after the elimination of the n shell. To first order of 
perturbation, we obtain 

{-iu + v^^)k^ + 6v^n)k^) uf{k) = -^Pi}"„(k) j dp[uf {p)u<{k - p) 

-b<{p)b<{k-p)] (11) 
(-zc^ + A(„)fc2 + a(„)A:2)6<(fc) = -iPr-^{k) J dp[u< {p)b^{k - p)] (12) 

where 

^-(«)(^) - (d- 1)P J,^,^, (2^)'^^*^^'^''^^.(„)(p)p2 + .(„)(,),2 

-Se{k,p,q) . . . ,2 ' I ^ 2 ] (13) 



(d-l)fc2 7p+^.fe (2^)'^^ ''^ '""'i^(n)(p)p2-|-A(„)(g)g2 

+S,{k,p,q) f f ^ 2 ] (14) 

with Si{k,p, q)s as functions of fc,p, and q. Hence, after the elimination of the {n + l)th shell, the effective viscosity 
and resistivity will be (i^, A)(„+i) (A;) = (t^, A)(„)(fc) + (5(t^, A)(„) (A;). 

We solve the above equations iteratively. To simplify, we substitute C{k) in the Eqs. (|l^, by one dimensional 
energy spectrum E(k) 

C^u,bb(^k) = (15) 

where Sj, is the surface area of d dimensional spheres. We assume that E'^{k) and E^{k) follow 

E'^ik) ^ i\:"H2/3fc-5/3; E^k) = E^'ikj/rA (16) 

Regarding and A(„) , we attempt the following form of solution 

(^,A)(„)(fc„fc') = iK-y/'u'/'k-*/'i^,X)l^ik') (17) 
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with k = kn+ik'{k' < 1) with the expectation that v*_^-^{k') and are universal functions for large n. We 

numerically solve for z/^*^^(fc') and A*^j(fc'). Our calculations reveal that the solutions oiv*^-^{k') and A*^^(fc') converge 

for all d> dc ~ 2.2. From this observation we can conclude that Kolmogorov's energy spectrum {E{k) oc k^^/^) is a 
consistent solution of the RG equations. Meanwhile, Kraichnan's 3/2 energy spectrum and vk'^ = Afc^ oc /cBq, where 
Bo (a constant) is the magnetic field of the large eddies, do not satisfy the renormalization group equations [Eqs. jl^ , 
p^)]. Hence E(k) oc k"^^'^ is not a consistent solution of RG equations. Our result regarding the nonexistence of 



stable RG fixed point for d = 2 is consistent with the RG calculation of Liang and Diamond |14|. Refer to Fig. (|l|) 
for illustration of and A^„-|(fc') for d = 3 and = 1. 

The values of renormalized parameters for d — 3 and various va are shown in Table |[ For large rA, the asymptotic 
h'* is close to the corresponding value for fluid turbulence, but the asymptotic A* is also comparable to v*. This implies 
that in fluid dominated regime, there is a signiflcant magnetic energy flux in addition to the usual Kolmogorov's flux 
in fluid modes. As is decreased, v* increases but A* decreases. This trend is seen till ta ~ 0.25, where the RG 
fixed point with nonzero v* and A* becomes unstable, and the trivial RG fixed point with i^* = A* = becomes 
stable. This result suggests an absence of turbulence for ta below 0.25. This is consistent with the fact that MHD 
equations become linear in the rA^O (fully magnetic) limit. 

We can proceed further and compute various cascade rates and Kolmogorov's constant for MHD using the renor- 
malized parameters computed above. To compute these quantities we resort to the energy equations, which are 

^ + '^'0 ^ (^-l)(2.W + kO X,,,,.o (^[^""("''P'^) + ^""("''^'P) 

+5«''(k'|p|q) + ^"^(k'|q|p)] (18) 

^ + ^ (d-l)(2.W + kO X,,,,.o (^[^"("''P'^) ^ ^"("''^'P) 

+5'"'(k'|p|q) + 5'"'(k'|q|p)] (19) 

where 

5™(k'|p|q) = -3([k'.u(q)][u(k').u(p)]), (20) 

5''^k'|p|q) = -S([k'.u(q)][b(k').b(p)]), (21) 

S-^(k'|p|q) = 3 ([k'.b(q)][u(k').b(p)]) , (22) 

5^"(k'|p|q) = -5«''(p|k'|q) (23) 

Here 3 stands for the imaginary part of the argument, and the above integrals have constraints that k' + p + q = 
(k = — k'). The energy equations in the above form have been written by Dar et al. [|6|, who interpret the terms 
S'(k|p|q) as energy transfer rate from mode p (second argument of S) to k (first argument of S) with mode q (third 
argument of S) acting as a mediator. This interpretation of energy transfer due to Dar et al. |]16[ is consistent with 
the earlier formalism. 

We can derive an expression for the energy transfer rate or energy flux from a wavcnumber sphere using S'(k'|p|q). 
The formula for the energy flux from inside of the X-sphere {X <) to outside of the F-sphere {Y >) is 

where X and Y stand for u or b. In our study we assume that the kinetic energy is forced at small wavenumbers, and 
the turbulence is steady. We calculate the above fluxes analytically to the leading order in perturbation series using 



the same procedure as that of Leslie 17 1. The flux is calculated using Eq. ( p4| ) by taking ensemble average of . 
The expression for (^S'''' {k\p\q)') is 



{S''''{k\p\q)) = I dt'[Ti{k,p,q)G''\k,t - t')C''\p,t - t')C''''{q,t ~ t') 
+Ts{k,p, q)G''\p, t - t')C''{k, t - t')C^^{q, t - t') 

+Twik,p, q)G^'\q, t - t')C'\k, t ~ t')C'\p, t - t')] (25) 
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where Ti{k,p,q) are functions of wavevectors k,p, and q. The expressions for other transfer rates {S'^'^{k\p\i 
(S'"''(fc|p|(7)), and (S'^"(A:|p|g)) look similar. In the above formulas we substitute Kolmogorov's spectrum [Eqs. 
for the energy spectrum, and the following expression for the effective viscosity and resistivity 

(i/,A)(/c) = (ii:")i/2ni/3fc-4/3(^*^^*) forfc>fc„ (26) 

Following the same procedure as Leslie p7|, we obtain the following nondimensional form of the equations 



^= {K^f^ f dv\n{l/v) I dw{vwf'^smaf-^F^< (27) 



n id-l)^Sd 







where the integrals Fy^ are function of u, w, v* , and A*. After a bit of manipulation we can obtain Hy^ /II and the 
constant . In addition we can also obtain the Kolmogrov's constant K for total energy 

E{k) = Kli^l^k'^/^ (28) 

using K = _K""(1 + r^). The values of Ily^/II and K for d = 3 and various ta are listed in Table |[ 

The entries in the Table | show that the cascade rates n|^> , n|j^, 11^^, 11^^ are approximately of the same order for 
rA between 0.5 and 1, but the flux IlJJ^is rather small. The sign of H^^is positive, indicating that the ME cascades 
forward, that is from large length-scales to small length-scales. The magnetic energy thus appearing at small length- 
scales will be lost due to resistive dissipation, and the large-scale magnetic field is maintained by the 11^^ flux. The 
Kolmogorov's constant K is approximately constant and is close to 1.6, same as that for fluid turbulence [ta = oo), 
all rA greater than 0.5. 

To summarize, we employed a self-consistent RG scheme to MHD turbulence and found that Kolmogorov's 5/3 
powerlaw is a consistent solution of RG equations for d > dc ~ 2.2. For Kolmogorov's solution, the renormalized 
viscosity and resistivity have been calculated, and they are found to be positive. For o? = 3, variation of v* and A* 
with TA shows some interesting features. For large r^, v* is same as that for fluid turbulence, but A* is also nonzero, 
in fact larger than v* . As ta is decreased, v* increases but A* decreases until rA ~ 0.25 at which value turbulence 
disappears. 

Using the flux equations we have obtained various fluxes and Kolmogorov's constant K. For d — 3, K does not 
vary significantly with the variation of r^, and it is close to K for fluid turbulence. We find that the cascade rate 
from magnetic-sphere to outside magnetic-sphere (H^^) is positive, a result consistent with the numerical results of 
Dar et al. |§. 

In this paper we have restricted ourselves to nonhelical turbulence. Helical MHD turbulence is very important 
specially in the light of enhancement of magnetic energy (dynamo). However, the physics of helical turbulence is more 
complex with the appearance of inverse cascade of magnetic helicity etc. The field-theoretic analysis for this case will 
be taken up later. Recent studies show that the mean magnetic field has a strong effect on energy spectrum, and it 
induces anisotropy. A full-fledge fleld theory calculation in the presence of mean magnetic field is also necessary for 
a clearer picture of MD turbulence 

The author thanks J. K. Bhattacharjee for very valuable discussions and ideas. He also thank G. Dar and V. 
Eswaran for many insights from numerical results. 



[1] R. H. Kraichnan, Phys. Fluids 8, 1385 (1965). 
[2] P. S. Iroshnikov, Sov. Astron. I. 7, 566 (1964). 

[3] E. Marsch, in Reviews m Modern Astronomy, edited by G. Klare (Springer- Verlog, Berlin, 1990), p. 43. 

[4] W. H. Matthaeus and Y. Zhou, Phys. Fluids B 1, 1929 (1989). 

[5] Y. Zhou and W. H. Matthaeus, J. Geophys. Res. 95, 10291 (1990). 

[6] M. K. Verma et al, J. Geophys. Res. 101, 21619 (1996). 

[7] W. C. Miiller and D. Biskamp, Phys. Rev. Lett. 84, 475 (2000). 

[8] D. Biskamp and W. C. Miiller, Phys. Plasma 7, 4889 (2000). 

[9] M. K. Verma, Phys. Plasma 6, 1455 (1999). 
[10] S. Sridhar and P. Goldreich, Astrophys. J. 432, 612 (1994). 
[11] P. Goldreich and S. Sridhar, Astrophys. J. 438, 763 (1995). 



4 



[12] D. Forstcr, D. R. Nelson, and M. J. Stephen, Phys. Rev. A 16, 732 (1977); C. DeDominicis and P. C. Martin, Phys. Rev. 
A 19, 419 (1979); J. D. Fournier and U. Frisch, Phys. Rev. A 17, 747 (1979); V. Yakhot and S. A. Orszag, J. Sci. Comput. 
1, 3 (1986). 

[13] W. D. McComb, The Physics of Fluid Turbulence (Claredon, Oxford, 1990); W. D. McComb, Rep. Prog. Phys. 58, 1117 
(1995). 

[14] J. D. Fournier, P.-L. Sulem, and A. Pouquet, J. Phys. A 15, 1393 (1982); S. J. Camargo and H. Tasso, Phys. Fluids B 4, 

1199 (1992); W. Z. Liang and P. H. Diamond, Phys. Fluids B 5, 63 (1993). 
[15] M. M. Stanisic, Mathematic Theory of Turbulence (Springer- Verlog, New York, 1988). 
[16] G. Dar, M. K. Verma, and V. Eswaran, Physica D, 157, 207, 2001. 
[17] D. C. Leslie, Development in the Theory of Turbulence (Claxedon, Oxford, 1973). 



TABLE I. The values of u*, A*, i^""*, i^"''*. A*""*, A*'*'* for for various rA when d = 3 and (Tc = 



rA 


z/* 


A* 


n:;$/n 


n;^</n 




n^$/n 


u-</n 


K 


oo 


0.38 




1 










1.53 


5000 


0.36 


0.85 


1 


3.5E-4 


-1.05E-4 


2.4E-4 


1.3E-4 


1.51 


5 


0.47 


0.82 


0.61 


0.26 


-0.050 


0.19 


0.13 


1.51 


1 


1.00 


0.69 


0.12 


0.39 


0.12 


0.37 


0.49 


1.50 


0.5 


2.11 


0.50 


0.037 


0.33 


0.33 


0.30 


0.63 


1.65 


0.3 


11.0 


0.14 


0.011 


0.36 


0.42 


0.21 


0.63 


3.26 


0.2 
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FIG. 1. Plot of v*{k') (solid lines) and A*(A;') (dashed lines) vs. k' for d = 3 and ctc = 0, = 1- Values at various iterations 
are shown by different curves. 
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